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Abstract 

A  non-Cartesian  coordinate  system  is  developed  which  permits  the  vibrational 
motions  of  Buckminsterfullerene  (Bucky  ball)  to  be  expressed  in  terms  of 
four  force  constants.  A  180  *  180  matrix  is  then  derived  which,  when 
diagonalized,  yields  the  complete  vibrational  spectrum.  These  results  are 
compared  with  those  obtained  previously  via  a  MNDO  calculation. 
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I.  Introduction 

The  purpose  of  this  paper  is  to  investigate  the  vibrational  motions  of 
Buckminsterfullerene  (Bucky  ball)  in  a  way  that  lends  itself  toward  future 
application.  Hence  simplicity,  both  theoretical  and  practical,  is  an 
important  consideration.  To  this  end  we  can  take  advantage  of  the  synanetry 
of  the  species,  and  we  shall  derive  a  relatively  simple  procedure  for 
calculating  the  vibrational  normal  frequencies.  Furthermore,  our  method 
lends  itself  to  physical  interpretation,  and  we  shall  be  able  to  make  some 
statements  about  the  nature  of  Bucky  ball.  The  experimental  literature  on 
Bucky  Ball  is  well  known. *  Theoretically  much  work  has  been  dedicated  to 
determining  the  most  stable  structure  of  clusters, ^  *®  and  to 
determining  some  of  their  optical  properties.**  We  continue  this  work  here 
by  assuming  that  Bucky  ball  exists  and  is  relatively  stable. 

In  essence  our  calculation  can  be  described  as  follows.  We  express  the 
position  of  each  atom  in  local  coordinates,  choosing  our  coordinate  system 
(non-Cartesian)  in  a  way  that  simplifies  the  calculation.  This  permits  us 
to  express  the  equation  of  motion  of  each  atom  in  terms  of  four  force 
constants,  and  since  all  atoms  are  identical  we  immediately  derive  an 


equation  of  motion  for  the  entire  system.  The  fundamental  principles  behind 


this  calculation  are  found  in  Ref.  12.  While  we  must  input  the  force 


constants  into  our  calculation  (this  is  a  disadvantage),  we  generate  the 


correct  number  of  modes,  the  correct  point  group  properties,  and,  at  least 


qualitatively,  the  correct  frequencies.  This  requires  only  the 


diagonal! sat ion  of  one  180  *  180  matrix,  a  procedure  which  requires  only 


about  three  minutes  CPU  on  a  VAX.  Furthermore,  since  we  can  vary  the  force 


13 

constants  and  since  we  can  compare  our  results  with  those  of  Stanton,  we 

are  able  to  make  some  predictions  about  the  structure  of  Bucky  ball  as 

compared  to  benzene.  In  the  next  section  we  derive  our  method  for 

calculating  the  normal  mode  frequencies,  and  in  Section  III  we  compare  our 

13 

results  with  those  of  Stanton  and  discuss  the  implications. 

II.  Theory 

The  structure  of  Buckminsterfullerene,  illustrated  in  Fig.  1,  is  that 

of  a  truncated  icosahedron  with  sixty  vertices,  twenty  hexagonal  faces  and 

twelve  pentagonal  faces.  A  carbon  atom  occupies  each  vertex.  The  bonds 

separating  a  hexagon  from  a  pentagon  are  found  to  be  more  "single"  than 

6  8 

bonds  separating  two  hexagons.  '  The  bond  lengths  are  given  as  1.41  A  and 

g 

1.54  A,  respectively.  Since  each  vertex  is  at  the  intersection  of  two 
hexagonal  and  one  pentagonal  face,  it  follows  that  it  is  joined  by  one 
"double"  bond  and  two  "single"  bonds.  This  implies  that  the  motion  of  an 


atom  around  it  equilibrium  position  can  be  decomposed  along  the  three  non- 

coplanar  directions  of  the  adjacent  bonds  (Fig.  2).  We  can  use  x^,  y^  and 

z^  to  denote  the  deviation  of  atom  i  from  its  equilibrium  position  in  these 

directions  (we  caution  the  reader  not  to  confuse  our  notation  with  Cartesian 

coordinates,  since  ours  are  not  Cartesian).  The  angle  between  two  adjacent 
3n 

single  bonds  is  j-  and  that  between  a  double  bond  and  an  adjacent  single 

bond  is  y4  One  can  work  out  the  angle  between  a  double  bond  and  its 

adjacent  pentagonal  face  as  T  ■  cos  *( - 4-).  The  angle  between  two 

Icosf- 

adjacent  hexagonal  faces  is  ♦  ■  cos  *(|(siny^)^  *  1],  and  the  angle  betwee 
adjacent  hexagonal  and  pentagonal  faces  is 
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The  variations  of  the  bond  length  a^,  and  (see  Fig.  2)  due  to 
the  motion  of  the  atoms  are,  to  first  order, 


4*ij  ’  (‘i  •  k  '  k>  +  (,j  '  k  '  kK 

4bik  '  -  2*1  -  cos5^* i ^  +  <‘k  •  K  •  co*r5'k)>  <l> 

,  /  1  2s  \  /  1  2s _  . 

6cu "  (*i  ■  2*i  ■  cosryi)  +  (yt  ■  2*t  ■  co*rV- 

Also,  the  variations  of  the  angles  between  bonds  a,  3  and  y  are 

6aj  »  £(-sin^  cos*  +  sin^  Zj  +■  sin^  x^  -  sinj^cos© 

-  sin^  y^  +  sin^  cos*  z^  -  sin^  +  sin|^  cos©  z^), 

6pj  •  £(-sin^  cos*  Zj  +  sin^  y^  +  sin^  x^  -  sin^cos©  z^ 

-  sin^  z A  +  sin^  cos*  y^  -  sin^  Xj  +  sin^  cos©  y^),  (2) 

6y^  ■  £(sin^  y^  -  sin^  cos©  x^  +  sin|^  z^  -  sin^  cos©  x^ 

♦  sin^  cos©  x^  -  sin|^  y^  +  sin^  cos©  x^  -  sinj^  z^), 

where  L  refers  to  the  bond  lengths  which  are  assumed  equal. 

The  kinetic  energy  of  each  atom  is 


(3) 


Ti  *  2t(xi  '  Fi  ‘  IV  +  <yicos5  '  *icos5) 
+(yicos|^  sinf  +  ^cosy^  sinf)2). 


where  o  is  the  mass  of  each  atom.  It  is  assumed  that  the  atoms  undergo 
harmonic  oscillations  around  their  equilibrium  positions,  where  the 
potential  is  due  to  the  variations  of  the  bond  lengths  and  the  angles 
between  them.  If  atoms  j,  k  and  t  represent  nearest  neighbors  of  atom  i 
as  shown  in  Fig.  2,  then  the  Lagrangian  for  this  system  takes  the  form 


L  '  T  ■  v  *  l  Ti  •  l  k4alj  -  l  2*2<bik  -  l 

i  i>j  i>k  i>t 


-  £  +  <*2)  *  VYih 


where  k^,  k2»  k^  and  k^  are  all  force  constants.  Equation  (4)  can  be 
rewritten  in  the  compact  form 


1DV 

1  m  2  \  (Tmn  X*  X"  ’  X*  x"). 


m,nsl 


and  the  appropriate  Euler-Lagrange  equation  is 


lOV 

X"  -  -  £  (T-1)"*1  VmlXl. 
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By  diagonalizing  the  180  *  180  matrix  T  *V,  we  get  the  squares  of  the 
eigenfrequencies  along  with  their  associated  eigenvectors. 

Several  cases  are  worth  discussing: 

(a)  kj  *  kj  ■  k^  ■  0.  Each  double  bond  with  force  constant  becomes 

1/2 

independent.  Then  one  obtains  only  one  non-zero  eigenfrequency,  (k^/m)  , 

with  a  30-fold  degeneracy. 

(b)  k^  ■  k^  *  k^  *  0.  The  system  reduces  to  12  independent  pentagons.  For 

each  pentagon  the  eigenfrequencies  are  (with  degeneracies)  0.89401  (1), 
1.00029  (2)  and  1.15166  (2),  with  All  other  eigenfrequencies  are 

zero. 

(c)  k3  ■  k^  *  0.  One  obtains  90  non-zero  eigenfrequencies  corresponding  to 
the  90  bonds.  Since  there  is  no  force  on  the  angles,  Bucky  ball  can  be 
squished  in  90  -  6  ■  84  different  ways.  The  remaining  6  degrees  of  freedom 
are  due  to  translation  and  rotation. 

(d)  In  the  general  case,  there  are  always  6  zero  eigenfrequencies 
corresponding  to  global  translational  and  rotational  motions.  There  are 
also  three  non-degenerate  modes. 

III.  Results  and  Discussion 

Table  I  contains  180  eigenvalues  representing  the  vibrational 

frequencies  of  Buckminsterfullerene.  The  values  of  k^  -  k^  are  those  that 

13 

best  seem  to  fit  the  data  of  Stanton,  which  was  calculated  by  the  MNDO 
method  with  the  gradients  evaluated  at  six  points  around  an  atom.  The  best 
fit  was  determined  by  matching  the  highest  and  lowest  frequencies^  and  also 
by  matching  the  three  non-degenerate  frequencies.  In  addition,  we  have 
insisted  that  the  ratio  between  k^  and  ^  b®  approximately  equal  to  the 


* 


>1 

a 


1 


i 


i 


i 


••  ■ . •. •.«  v^j 
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inverse  ratio  of  the  bond  lengths.  Since  at  this  point  we  have  no 
independent  way  of  determining  k^,  we  have  insisted  that  it  equal  k^. 

For  benzene  the  appropriate  values  would  be  k^  ■  k^  ■  7.62  ■  10^ 
dynes/cm  and  k^  ■  .667  *  10^  dynes / cm.  ^  Comparing  these  values  with  those 
given  in  Table  I,  we  see  that  the  force  constants  for  Bucky  ball  are 
considerably  higher  than  for  benzene.  This  is  physically  reasonable  since 
the  stretching  of  a  bond  in  Bucky  ball  involves  a  dislocation  of  the  entire 
molecule,  and  hence  greater  resistance.  Similarly  we  expect  the  stretching 
constant  to  be  larger. 


It  is  also  instructive  to  look  at  the  three  non-degenerate  modes. 
Stanton^  has  calculated  their  frequencies  (in  cm  *)  as  611  (for  the  A 

8 

mode),  973  (Ay)  and  1667  (A  ).  Our  calculation  reveals  frequencies  at  548, 

970,  and  1627  cm  *,  respectively.  We  have  found  that  the  two  A  modes 

8 

depend  only  on  constants  k1  and  kj.  The  first  is  a  "breathing"  motion 
whereas  the  second  involves  oscillation  between  the  "single"  and  "double" 


bonds.  In  both  cases  angles  can  be  expected  to  play  no  role.  The  A^ 
motion,  on  the  other  hand,  involves  a  twisting  (angle -dependent)  motion  of 


the  pentagons.  When  k^  and  k^  are  zero,  the  eigenfrequency  of  this  motion 

becomes  zero.  All  other  modes  must  be  degenerate  according  to  the  point 

13 

group  of  Buckminsterfullerene. 


Our  method  is  well  suited  to  elucidating  this  sort  of  physical 
description  of  the  motions  since  the  k's  have  obvious  physical  significance. 
This  is  a  big  advantage  of  this  calculation,  along  with  the  numerical 
simplicity.  We  intend  to  use  this  method  to  further  investigate  the 
vibrational  spectrum  of  Buckminsterfullerene. 
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Frequency 

Symmetry 

Frequency 

Symmetry 

0.0 

- 

272.0 

Hg 

354.7 

Hu 

361.8 

F2u 

373.7 

Gu 

427.9 

Hg 

455.8 

Gg 

491.3 

Flu 

491.6 

Hu 

525.2 

F2g 

530.2 

Gg 

547.6 

Ag 

550.9 

Flu 

552.4 

Hg 

566.8 

Fig 

577.7 

Hu 

626.6 

F2u 

672.7 

Gg 

701.9 

Gu 

726.4 

F2g 

755.8 

Gu 

770.3 

Hu 

779.8 

Hg 

810.9 

Fig 

926.6 

F2g 

958.3 

Gu 

1019.0 

F2u 

1084.3 

Au 

1160.0 

Hg 

1173.5 

Gg 

1289.7 

Hu 

1309.4 

F2u 

1374.1 

Flu 

1398.5 

Hg 

1463.8 

F2g 

1578.5 

Hu 

1590.0 

Gg 

1620.4 

Gu 

1627.4 

Ag 

1655.3 

Flu 

1665.1 

F2g 

1688.2 

Hg 

1720.1 

F2u 

1764.4 

Gu 

1765.4 

Gg 

1830.0 

Hu 

1830.7 

Hg 

Table  1.  Frequencies  in  cm  *  and  their  degeneracies  for  the 
vibrational  modes  of  Buckminsterfullerene .  *  1.1  *  10^, 

kj  *  1.0  *  10^,  *  1.0  x  10^  and  k^  »  1.0  x  10^  dynes/cm. 


v'.v.  /  ■  .•  ^ cLUi  »v» ■ 
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Figure  1. 


Figure  2. 


FIGURE  CAPTIONS 

The  CgQ  cluster,  better  known  as  Buckminsterfullerene, 
is  a  truncated  icosahedron  with  20  hexagonal  faces  and 
12  pentagonal  faces. 

A  diagram  of  the  local  coordinate  system  used  in  the 
calculation. 
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